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ABSTRACT 

DEM analysis is a major diagnostic tool for stellar atmospheres. But both its derivation and 
its interpretation are notably difficult because of random and systematic errors, and the inverse 
nature of the problem. We use simulations with simple thermal distributions to investigate the 
inversion properties of SDO/AIA observations of the solar corona. This allows a systematic 
exploration of the parameter space and using a statistical approach, the respective probabilities 
of all the DEMs compatible with the uncertainties can be computed. Following this methodology, 
several important properties of the DEM inversion, including new limitations, can be derived and 
presented in a very synthetic fashion. 

In this first paper, we describe the formalism and we focus on isothermal plasmas, as building 
blocks to understand the more complex DEMs studied in the second paper. The behavior of the 
inversion of AIA data being thus quantified, and we provide new tools to properly interpret the 
DEM. We quantify the improvement of the isothermal inversion with 6 AIA bands compared to 
previous EUV imagers. The maximum temperature resolution of AIA is found to be 0.03 logTe, 
and we derive a rigorous test to quantify the compatibility of observations with the isothermal hy- 
pothesis. However we demonstrate limitations in the ability of AIA alone to distinguish different 
physical conditions. 

Subject headings: Sun: corona - Sun: UV radiation 



1. Motivation 

The Differential Emission Measure (DEM) di- 
agnostic technique offers crucial information about 
the thermal structuring of the solar and stellar 
atmospheres, providing a measure of the temper- 



ature distribution of the plasma along the line of 
sight (LOS). However, to derive the DEM from 
a set of observations is a complex task, due to 
the inverse nature of the problem, and the un- 
derstanding of its robustness and accuracy is still 
relevant today (e.g. Landi et al. 12011 : Testa et al 
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2012I ). Spectrometers are by nature better suited 
to DEM analysis than broad band imagers. But, 
because these latter generally offer a higher signal 
to noise ratio over a larger field of view (FOV), 
DEM codes have nevertheless been applied to the 
three coronal bands of the Extreme-ultraviolet 



Imag ing Telescope (EIT) (|Delaboudiniere et al. 
19951) or the Tra nsition Region and C oronal Ex 



plorer (TRACE, iHandv et al~lll999l ). However, 
these instruments were shown not to constrain 
the DEM enough to reach conclusive results. In 
recent years, the multiplication of passbands in in- 
strumen ts such as the X- Ra y Telescope (XRT) on 
Hinode ( Golub et al. 1120071 ) and th e Atmospheric 



Imaging Assembly (AIA) telescope (jLemen et al. 



20121 ) has brought new prospects to reliably es- 
timate the DEM simultaneously over a large 
FOV. Case studies of the properties of the in- 
version using the se instruments have be e n pub- 



lis hed by 
or 



Reale et al. 



e.g. , iMartinez-Svkora et al. I (|201l[ ) 



(120091) . 

Building on these results, the central objective 
of the work presented in this series of papers is to 
provide a systematic characterization of the DEM 
reconstruction problem to assess both its accuracy 
and robustness. Using our technique, the capabil- 
ities of a given instrument can be evaluated, and 
new tools facilitating the DEM interpretation are 
presented. We illustrate our methodology in the 
specific case of the six coronal bands of AIA, but 
the same principle can be applied to any set of 
broad band or spectroscopic measurements. 

Initially introduced for element abundance 



measurements, then furthe r develop ( 
(|l972l ) and I Jordan I 



Jefferies et al. 



ed by e.g., 
(|l97fih . the 

DEM formalism has been extensively used in the 
past decades, on most types of coronal structures , 



such as polar co ronal holes d Hahn et alj 120111 ) , 



polar plu mes (e.g. Del Zann a et al...2005l) ~stream- 



ers (e.g. Parenti et al.l 2000l). prom inences (e.g. 
Wiik et a l."l993':'Parenti fc Vial"2007l) quiet sun fe.g. 
Landi fc Landini 199& Parenti fc Viall HoOJ), 



bright poin ts (iBrosius et al 



gions (e.g. IWarren et al. 



20081) 
201lh . 



or active re- 
The thermal 



structuring of the stellar coronae is also investi- 



gated using the DEM analysis (e.g. lSanz-Forcada et al. 



20031) . In particular, the DEM is one of the tools 
commonly used to study the thermal stability of 
the coronal structures just mentioned, and to di- 
agnose the energy source balancing the observed 



radiative losses. For example, it can help to dis- 
criminate between steady or impulsive heating 
models predi cting different loop thermal struc- 
tures (see e.g . Klimchuk 2006 : Reale et al. 20101 : 



Susino et al.ll2010l: IWinebarger et al.l 120111 ). One 



of the approaches is to establish the cross field 
thermal structure of resolved loops which is then 
compared to the DEM simulated for unresolved 
multi-stranded and monolitic loops, impulsively 
or steadily heated. 

But reliably inferring the DEM from obser- 
vations has proved to be a genuine challenge. 
The fundamental limitations in the DEM inver- 
sion h av e been discussed b y, e.g. . 



Jefferies et al. 



19721) : ICraig fc Brown I ([19761) : iBrown et al 



199l[ ): lJudge et all (|l997t ). including measure- 
ment noises, systematic errors, the width and 
shape of the contribution functions, and the as- 
sociated consequences of multiple solutions and 
limited temperature resolution. Many DEM inver- 
sion algorithms have been proposed to cope with 
these limitation s, each with its own strengths and 



weaknes ses (e.g. iWith broe 1975^; Craig fc Brown 



198 6: Judge et al. Landi fc Landini 199 
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K ashvap fc Drake lll99ilMcIntosh ll2000HWeber et al 



2004 iGorvaev et al. I l2010t iHannah fc Kontar 
2012| ). In parallel to these developments, au- 
thors have been attentive early on to estimat e 
the accuracy of the inversions (e.g. [Dere][i97i), 
eventually comparing se veral algorithms (e.g. 
Fludra fc Svlwester Ill986l ). 

Due to the intrinsic underconstraint of inverse 
problems and to the inevitable presence of ran- 
dom and systematic measurement errors, multiple 
physical solutions consistent with the observations 
exist, even if mathematical uniqueness and sta- 
bility can be ensured via, e.g., regularization. It 
is nevertheless possible to quantify the amount of 
knowledge, or ignorance, on the physical parame- 
ter of interest by rigorously defining levels of con- 
fidence in the possible solutions or classes of solu- 
tions that can explain the observations within the 
uncertainties. This is a desirable feature for any 
inversion scheme if it is to be able, for example, to 
discriminate or even to define, isothermality and 
multithermality. 

In this perspective, we developed a technique 
to systematically explore the whole space of solu- 
tions, in order to determine their respective prob- 
abilities and quantify the robustness of the inver- 
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sion with respect to plasma parameters, random 
and systematic errors. We used data simulated 
with simple DEM forms to systematically scan a 
wide range of plasma conditions, from isothermal 
to broadly multithermal, and several inversion hy- 
potheses. Comparing the DEM solutions to the 
input of the simulations, it is possible to quan- 
tify the quality of the inversion. Following this 
strategy, we are able to completely characterize 
the statistical properties of the inversion for sev- 
eral parametric DEM distributions. We argue that 
even though the specifics may vary, the main con- 
clusions concerning the existence of multiple solu- 
tions and the ability to distinguish isothermality 
from multithermality also apply to more generic 
forms of DEM distributions. 

In this first paper, we focus on the response of 
AIA to isothermal plasmas. The properties of the 
isothermal inversion thus observed will serve as 
building blocks for the interpretation of the more 
complex DEM solutions studied in the second pa- 
per (hereafter Paper II) . Section [2] describes the 
general methodology and the practical implemen- 
tation in the case of AIA, including the data simu- 
lation, the inversion scheme, the sources of random 
and systematic errors, and the different DEM dis- 
tribution models considered. Results for isother- 
mal plasmas are presented and discussed in Sec- 
tion [3l A summary introducing the treatment of 
more generic DEM forms is given in conclusion. 

2. Methodology 
2.1. DEM formalism 

Under the assumption that the observed plasma 
is optically thin, integration along the line of sight 
(LOS) of collisional emission lines and continua 
produces in the spectral band b of an instrument 
an intensity 



h = 



1 

An 



Rb{ne,Te)nlds (1) 



where Rb{ne,Te), the response of the instrument 
to a unit volume of plasma of electron number 
density rig and temperature Tg, is given by 



Rhine, Te) ^^SbiXl)Ax Gxj{^e,Te) 

x.i 
Jo 



(2) 



The first term of the right member accounts for 
each spectral line / of each ionic species X of abun- 
dance Ax, and the second term represents the 
contribution of the continua. 5*6 (A) is the spec- 
tral sensitivity of the band b of the instrument. 
The respective contribution functions Gx,i{ne, Te) 
and Gc{ne,Te) of the lines and continua con- 
tain the p hysics of the radiation emission pro- 
cesses (e.g. iMason fc Monsignori Fossi II1994I ) and 
can be computed using the relevant atomic data. 
As long as one considers total line intensities, 
equations ([T]) and ([2]) are generic and apply to 
imaging telescopes as well as to spectrometers. 



Summarizing the original reasoning of lPottasch 

(|l963 , 1964). since the function Rf,{ne,Te) is gen- 
erally weakly dependent on the density and is 
peaked with temperature, gives a measure of 
J n^ds where the integration is now limited to 
tfie portions p of the LOS where the temperature 
is such that significant emission is produced. If 
measurements are available at several wavebands, 
it is possible to plot nids as a function of the 
bands' peak temperatures. Generalizing this logic 
into a differential form, and assuming that the el- 
ement abundances are constant, equation ([T|) can 
be reformulated as 



1 

47r 



Rb{Te)i{Te)d\0gTe 



(3) 



where i(Te) = nl (Te) dp /d log Te is the DEM, that 
provides a measure of the amount of emitting 
plasma as a function of tem p eratu r j^. As demon- 
strated bv ICraig fc Brown I (|l976l) . nf(Te) is the 
mean square electron density over the regions dp 
of the LOS at temperature Te, weighted by the in- 
verse of the temperature gradients in these regions. 
The total Emission Measure (EM) is obtained by 
integrating the DEM over the temperature 



EM = I ^{Te)d\ogTe= / nlds. 
Jo Jo 



(4) 



Solving the DEM integral equation ([3|) implies re- 
versing the image acquisition, LOS integration and 



^The logarithmic scale is justified by the shape of the con- 
tribution functions (see Figure [8ll. The DEM can also be 
defined in linear scale as 5(Te) = n'^{Te)dp/dTe. There is a 
factor dlogTe/dTe = 1/lnlO Tg between the two conven- 
tions. 
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photon emission processes to derive the distribu- 
tion of temperature in the solar corona from ob- 
served spectral line intensities. We will now inves- 
tigate the properties of this inversion. 

2.2. Probabilistic interpretation of the 
DEM solutions 

Let us consider a plasma characterized by a 
DEM ^^(Te). The corresponding intensities ob- 
served in Nb spectral bands are noted Ih^^i^,^)- In 
order to solve the DEM inverse problem - estimat- 
ing from the observations - one uses a criterion 
C(^) that defines the distance between the data 
I^''^ and the theoretical intensities /^^(C) com- 
puted using equations ([2]) and ([3]) for any DEM 
^(Te). By definition the DEM ^^(re) solution of 
the inversion is the one that minimizes this crite- 



rion: 



= argminC(^). 



(5) 



Since the 1^^^ are affected by measurement noises 
and the I^'^ by systematic errors in the calibra- 
tion and atomic physics, the inversion can yield 
different solutions of probabilities P{£,^\£,^) for 
a given DEM of the plasma. Bayes' theorem 
then gives 



') 



(6) 



which is the conditional probability that the 
plasma has a DEM knowing the result of 
the inversion. P(^^) = / P{C^^^)P{^^) d^^ is 
the total probability of obtaining whatever . 
In the Bayesian framework, P{£,^) is called the 
prior. It is uniformly distributed if there is no a 
priori information on the DEM of the plasma. 
Conversely, a priori knowledge or assumptions on 
the plasma are represented by a varying P(^^). 
For example, zero probabilities can be assigned to 
non physical solutions. 

P{^^\£,^) contains all the information that can 
be obtained from a given set of measurements on 
the real DEM of the plasma and as such, it 
is a desirable quantity to evaluate. Indeed, if the 
DEM is to be used to discriminate between physi- 
cal models, as it is for example the case in the coro- 
nal heating debate, finding a solution that mini- 
mizes the criterion is necessary, but it is not suffi- 
cient. It is also crucial to be able to determine if 



other solutions are consistent with the uncertain- 
ties, what are their respective probabilities, and 
how much they differ from each other. 

In principle, and without a priori on the 
plasma, P{^^\£,^) and thus P{^^\£,^) can be esti- 
mated for any minim ization scheme using Monte- 
Carlo simulations (jMetropolis fc Ulam I 



For each the Nb observed ''(C^) are sim- 
ulated using equations © and ^ and adding 
photon and instrumental noises. Systematic er- 
rors are incorporated to the and the resulting 
criterion is minimized. P{£,^\^^) is then evalu- 
ated from the N solutions corresponding to N 
realizations of the random variables. But since 
several can potentially yield the same the 
derivation of P{^^\£,^) from equation ^ requires 
to know P{£,^), the probability to obtain what- 
ever This is generally not possible, for it 
requires the exploration of an infinite number of 
plasma DEMs. 

This is why DEM inversion research often fo- 
cuses on the minimization part of the problem, 
P{i^\£,^) being supposed to be well behaved be- 
cause of the proper choice of prior and the multi- 
plication of passbands or spectral lines. However, 
P{^^\(^) can be computed if the DEM of the 
plasma can be described by a limited number of 
parameters. In this case, one can scan the whole 
parameter space and use the Monte-Carlo simula- 
tions to estimate P{£,'\^^) for all possible . The 
possibility that multiple yield an identical in- 
version solution being now taken into account, 
one can determine P{£,^) and thus derive P{£,^\(,^) 
from equation ([5]). 

This limitation of the complexity of the DEMs 
that can be considered corresponds to adopt- 
ing a non- uniform prior P{£,^), while proba- 
bilistic treatments were justly developed with 
the opposite objective of relaxing such non- 
ph ysical assumptions (e.g. the MCMC method 
of iKashvap fc Drake mm . But rather than the 
development of a generic DEM inversion method, 
our objective is to study the behaviour of P(^^|^^) 
in controlled experiments. And if the parameteri- 
zation is properly chosen, the can still represent 
a variety of plasma conditions, from isothermal to 
broadly multithermal. In addition, we did not 
make any assumption on the number and prop- 
erties of the spectral bands, nor on the defini- 
tion of the criterion nor on the algorithm used to 
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minimize it. The method described to compute 
P{£,^\£,^) can therefore be used to characterize any 
inversion scheme in the range of physical condi- 
tions covered by the chosen distributions. 

2.3. Inversion method 

Devising an efficient way to locate the abso- 
lute minimum of the criterion is not trivial. For 
example, without further assumption, its defini- 
tion alone does not guarantee that it has a single 
minimum, so that iterative algorithms may con- 
verge to different local minima depending on the 
initial guess solution. Furthermore, if the value of 
the minimum itself is a measure of the goodness 
of fit, it does not provide information on the ro- 
bustness of the solution. How well the solution is 
constrained is instead related to the topography of 
the minimum and its surroundings; the minimum 
may be deep or shallow and wide or narrow with 
respect to the different parameters describing the 
DEM curve. 

The number of DEMs resulting in significantly 
different sets of intensities within the dynamic 
range of an instrument is potentially extremely 
large. However, a systematic mapping of the cri- 
terion aimed at revealing its minima and their to- 
pography is possible if the search is restricted to 
a subclass of all possible DEM forms. Indeed, if 
the DEM is fully determined by a limited number 
of parameters, one can regularly sample the pa- 
rameter space and compute once and for all the 
corresponding theoretical intensities The 
criterion, i.e. the distance between the II^ and 
the measured J^''^, is thus computable as a func- 
tion of the DEM parameters for any given set of 
observations. It is then trivial to find its absolute 
minimum and the corresponding DEM solution 
or to visualize it as a function of the DEM param- 
eters. 

2.4. Implementation 

The procedure used to compute P(C'^IC^) is 
summarized in Figure [TJ The parametric DEM 
forms are described in section [2.4. II The intensi- 
ties I^^" observed in bands are the sum of av- 
erage intensities 1° and random perturbations Uh 
due to photon shot noise and measurement errors 

it^lUnb. (7) 



The are equal to the theoretical intensities /^^ 
in the case of a hypothetically perfect knowledge 
of the instrument calibration and atomic physics. 
In practice however, the if'' are affected by sys- 
tematic errors Sb 

It = 4° + Sb. (8) 

Since there is no way of knowing whether the in- 
tensities that can be computed from equations ([2]) 
and ([3]) for any DEM ^ are overestimated or un- 
derestimated, we identify thenjl to the reference 
theoretical intensities Ijj . The distributions of ran- 
dom and systematic errors are discussed in sec- 
tion 12.4.31 The detail of the calculation of the 
/° is given in section 12.4.21 From these, we can 
either simulate observations I^^^ by adding mea- 
surement noises Ub (equation (O), or obtain var- 
ious estimates of the if by adding perturbations 
representing the systematics Sb (equation ([5])). 

The criterion C(^) and the corresponding mini- 
mization scheme are described in section [2.4.4l For 
any plasma DEM <^^, Monte-Carlo realizations of 
the noises Ub and systematics Sb yield several es- 
timates from which we compute P{£,^\£,^)- Fi- 
nally, P{£,^\^^) is obtained after scanning all pos- 
sible plasma DEMs (section [2A5]). 

2.4-1- DEM distribution models 

Ensuing the discussions of sections 12.21 and 12.31 
the and are both constrained to belong to 
one of the three following classes of DEM distri- 
butions defined by two or three parameters: 

• Isothermal 

^,,oiTe)^ EM SiT,-T,), (9) 

where the DEM is reduced to a Dirac S func- 
tion centred on the temperature Tc- EM is 
the total emission measure defined by equa- 
tion (HI). 

• Gaussian in log 

igau{Te) =EMN{\ogT, - logT,), 

with TV (a;) ^ ^ cxp -— 7 

^It is also possible to adopt the view that the intensities 
computed with CHIANTI are one of the possible estimates 
of the if, in which case we obtain the /j^ by adding sys- 
tematic errors. The only difference between the two con- 
ventions is the sign of Sf,. The criterion and therefore the 
results are identical in both cases. 
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The plasma is here predominantly dis- 
tributed around a central temperature Tc 
with a width a. 



• Top hat in log Tg 

aat(Te) ^EM n(l0gTe 

if \x\ < 



with n(x) = 



(11) 



else 



The plasma is uniformly distributed over a 
width a around Tc- 

There is no reason for the solar plasma to follow 
one of these distributions, nor are they the only 
possible choices. But even though they are simple 
enough to allow a detailed analysis of the prop- 
erties of the DEM inversion, they can nonetheless 
represent a variety of plasma conditions. The con- 
clusions drawn can therefore help understand the 
behaviour of more generic DEM forms. Further- 
more, since the class of solution DEMs does not 
have to be the same as that of the plasma DEMs 
it is possible to investigate the impact of a 
wrong assumption on the shape of the DEM. For 
example, one can compute for isothermal 

solutions while the plasma DEM is multi- 
thermal (see paper II). 

2.J^.2. Reference theoretical intensities 

Equations © and ^ are used to compute the 
reference theoretical intensities /°(^) for any DEM 
5. They are then used to form both simulated ob- 
servations and various estimates of the theoretical 
intensities with equations ^ and 

From equations (0), pUj) and (ITT]) , we derive 
the expressions of these reference intensities as a 
function of the parameters EM, Tc and a for the 
three types of DEM distributions. 



Isothermal 

I^{EM, Tc) = EM / Rb{Tc) S{Te ~ d 
Jo 

= EMRb{Tc) 

(12) 

Gaussian 

r+00 

I^iEM, Tc, a) = EM / Ri,{Tc) AA(log - 
Jo 

= EM{Rb*J^){Tc,a) 

(13) 



logTe 



• Top hat 

I°{EM,Tc,a) 



= EM / RhiTc) n(log Tc - log Tc) d log Tc 
Jo 

= EM{Rb*U){Tc,rj) 



(14) 



We note that in all cases, the reference theoreti- 
cal intensities are equal to the convolution prod- 
uct of the instrument response function Rb{Te) by 
the chosen DEM S,{Te). The /° are pre-computed 
for all possible combinations of parameters EM , 
Tc, and cr. The appropriate range and resolution 
to be used for each parameter can be determined 
from plausible plasma properties and taking into 
account the instrument characteristics. 

The responses Rb{Te) of the six AIA coronal 
bands are computed using equation The 
contribution functions G{Tc) are obtained us- 
ing the v ersion 7.0 of the CHIA NTI atomic 
database (|Dere et al. I Il997l l2009l) . We used 
the CHIANTI ionization balance and the ex- 
tended coronal abundances. The summation 
is extended over the 5 nm to 50 nm spectral 
range for all bands. The instrument sensitiv- 
ity Sb{X) is obtained as a function of wave- 
length in units of DN.cm^.ph~"'^.sr~^ by calling 
the function aia_get jresponse provided in the 
AIA branch of the Interactive Data Language 
(IDL) Solar Software (SSW) package with the /DN, 
/area and /full keywords. This function imple- 
m ents the AIA pre-fligh t calibration as described 
in lBoerner et al. 1 ( 2012h . Since photon shot noise 
must be taken into account in the error budget 
(section 12.4.31) , the /° (^) must be computed for 
given exposure times and not per second. We 
used the standard AIA exposures of 2 s for the 
17.1 nm and 19.3 nm bands, and 2.9 s for the 
others. 

The contribution functions are computed us- 
ing CHIANTI from log(re) = 5 to log(Te) = 7.5 
in steps of 0.005 log(Te), oversampling the CHI- 
ANTI grid by a factor 10 using cubic spline in- 
terpolations. The emission measure varies over 
a wide range from 10^^ cm~^ to 10^'^ cm~^ in 
steps of 0.04 log(£'M). The DEM width varies lin- 
log Tc)^Fi^ i'l steps from cr = to cr = 0.8 log(re). 
This choice of sampling leads to pre-computing 
10^ groups of 6 AIA intensities, which represents 
easily manageable data cubes. 
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2.4-3. Uncertainties 

Uncertainties due to random and systematic er- 
rors are at the heart of the problem of the DEM 
inversion. The two affect the observations and 
their interpre tation in different manners (see e.g. 



Taylor 1 119971 ). Observations are mostly affected 



by random errors caused by both Poisson photon 
shot noise and nearly Gaussian detection noises 
like thermal and read noise. These noises vary 
randomly from pixel to pixel and from exposure to 
exposure. On the other hand, the errors made on 
the calibration and atomic physics systematically 
skew the interpretation of all observed intensities 
by the same amount and in the same direction. 

It is possible to realistically simulate in the I^''* 
the statistical properties of the noises affecting the 
data. The reference intensities have units of 
Digital Numbers (DN). The number of electrons 
collected in each pixel over the exposure time is 
obtained by multiplying these values by the gains 
(in e~/DN) of the detectors' analog to digital con- 
verters listed in SSW. The number of detected 
photons is then obtained by dividing the result 
by the quantum yield of the detector, i.e. the 
number of photoelectrons produced per interact- 
ing photorlfl. These photon intensities are then 
perturbed by Poisson noise and converted back to 
photo electrons. 22 e " RMS of Gaussian CCD read 



noise (jBoerner et al. 1120121 ) are finally added be 



fore conversion to DN. 

Determining the statistical properties of the 
systematic errors is more challenging. The tab- 
ulated calibration and atomic physics provides a 
single estimate of the instrument response Rb, but 
systematics nonetheless have a probability distri- 
bution. Indeed, the calibration is the result of 
laboratory measurements themselves affected by 
random and systematic errors. If we could recali- 
brate the instrument a number of times in different 
facilities we would obtain a distribution of instru- 
mental sensitivities Sb{X), the adopted calibration 
corresponding to one of them. Likewise, different 



^An approximation of the quantum yield of silicon is given 
by he/ {3.65 q A) where 3.65 is the energy in eV required to 
create an electron hole pair, q is the elementary charge, c 
is the speed of light in vacuum and h is Planck's constant. 
Note that in this calculation we assume that all interacting 
photons have the same wavelength. However, since the 
full width at half maximum of the AIA bands is comprised 
between 0.2 and 1.0 nm, the error made is only a few 10^'^. 



atomic physics codes will give different estimates 
of the contribution functions G(ne,Te), the CHI- 
ANTI output being one of them. It is however 
difhcult to characterize these two probability dis- 
tributions. They are generally implicitly assumed 
to be Gaussian and the adopted values to be the 
most probable. But the distributions may in fact 
be uniform, or asymmetric, or biased, etc. 

The calibration involves a complex chain of 
measurements, the uncertainties of which are dif- 
ficult to track and estimate. After independent 
radiometric calibrations, comparable EUV instru- 
ments on S OHO were found to agree only within 
about 25% (iHuber et al. Il2002h . Subsequent com- 
parisons could not resolve the discrepancies nor 
identify their origin in random errors or biases in 
the individual calibrations. We can only say that 
the adopted calibration of every SOHO instrument 
introduces a systematic error in the data analysis 
but without being able to tell how much and in 
what direction. It is likely that inter-calibration 
between AIA and other instruments would run 
into similar limitations. 

Errors in the contribution functions are a major 



contr i butor to the uncert ainties (e.g. iLang et al 



Judge et al. I 119971 ) 



Since the proper- 
transitions are de- 



1990t 

ties of the known atomic 
rived either from measurements or modelling, 
they are not infinitely accurate. Missing tran- 
sitions lead to underestimated contributions func- 
tions, as it is t he case for the 9.4 nm chan- 
nel of AIA (e.g. Aschwanden fc B ocrner ' 12011 : 
O'Dwver et al. ll201lHFoster fc Testa . 2011,) . The 



abun dances are affected by about 10% uncertain- 
ties ( Asplund et al. 2009() . not taking into ac- 
count p ossible local e nhancements of high FIP el- 
ements (| Young 1120051 ). These imply that, at least 
in some cases, the abundances are not constant 
along the line of sight, as assumed in the DEM 
analysis. The plasma may not be in ionization 
balance, in which case the CHIANTI calculations 
of transition rates are not valid. The response 
functions i?f, of AIA are also not independent 
from the electron number density, which is one 
of the assumptions made in deriving the DEM 
expression from equations ([T|) to ([3]). When us- 
ing spectrometers, the spectral lines are chosen so 
that this hypothesis is effectively verified. We plot 
in Figure[5]the normalized maximum of Ri,{T(,, n^) 
versus electron number density. In the AIA field of 
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view, rie can vary fro m about 10^ cm ^ in corona l 
holes at 1.2 R© (e.g. 



Guhathakurta at al. I[l999h 



to about 10 ^° cm ^ in dense coronal loops (e.g. 
lRealell2002l) . In this range, only the 9.4 nm band 
(solid line) is completely independent on the den- 
sity. The response function of all other bands 
decreases as the density increases, the variation 
reaching about 35% for the 17.1 nm band (short 
dashed line). Since the contribution functions 
have to be computed for a constant electron num- 
ber density (we chose 10^ cm~'^), they are re- 
spectively under or over-estimated if the observed 
structures are more or less dense. The impact 
can be mitigated if one has independent knowl- 
edge of the range of densities on the LOS, but 
it nonetheless represents an additional source of 
uncertainty compared to using density insensitive 
spectral lines. Finally, these various sources of 
uncertainties do not affect all spectral bands by 
the same amount. 

Rigorously estimating the properties of the 
probability distributions of the systematic er- 
rors would thus require a detailed analysis of 
the calibration process and of the atomic physics 
data and models that is beyond the scope of 
this paper. In these conditions, we make the 
simplifying assumption that all systematics are 
Gaussian distri b uted and unbiased. According to 
iBoerner et al. ( 2012 ). uncertainties on the pre- 
flight instrument calibration are of the order of 
25%. This is thus interpreted as a Gaussian prob- 
ability distribution centred on the tabulated val- 
ues with a 25% standard deviation. Likewise, we 
used 25% uncertainty on the atomic physics for 
all bands, typical of the estimates found in the 
literature. Calibration and atomic physics uncer- 
tainties were added quadratically for a net 35% 
uncertainty on the response functions Rb. The 
I^'^ are thus obtained by adding Gaussian random 
perturbations to the 

2.4-4- Criterion and minimization 

Since instrumental noises and systematic errors 
are assumed to be Gaussian distributed, we use a 
least square criterion 



E 

6=1 



(15) 



normalized to the total standard deviation of 
the uncertainties in each band, cr^ is obtained by 
summing quadratically the standard deviations of 
the four individual contributions: photon noise, 
read noise, calibration and atomic physics (sec- 
tion I2.4.3P . The value of the minimum of C (^) 
corresponding to the solution is noted 



= min C(0. 
From equations ^ and ^ we get 



E 

6=1 



(16) 



(17) 



If the family of solutions (Dirac, Gaussian or top 
hat) is identical to that of the plasma DEM 
then in the absence of noise = ^-nd the solu- 
tion given by equation ([5]) is strictly equal to 
However, in the presence of random and sys- 
tematic errors or if the assumed DEM form differs 
from that of the observed plasma, is not likely 
to be zero and the corresponding may be differ- 
ent from , for random fluctuations of Ub and Sb 
can compensate a difference between /°(C^) and 
IbiO- discussed in section 1531 properly inter- 
preting the value of provides a means of testing 
the pertinence of a given DEM model. 

Folding equation ^2^, or into equa- 
tion ([T5|) . we obtain the expression of C(^) for the 
corresponding DEM distributions. Given a set of 
observed intensities and a DEM model, the crite- 
rion can therefore be easily computed for all possi- 
ble combinations of the parameters i?M, Tc, and 
a using the /^(C) tabulated as described in sec- 
tion l2.4.^ Finding its minimum and thus the solu- 
tion is simplified to the location of the minimum 
of the C(^) matrix. This minimization scheme is 
not fast compared to, e.g., iterative gradient al- 
gorithms, but it ensures that the absolute mini- 
mum of the criterion is found whatever its topog- 
raphy. Furthermore, this operation can be effi- 
ciently implemented on the Graphics Processing 
Units (GPU) of modern graphics cards by using 
their CUDA capability. We implemented a scheme 
in which each GPU core is in charge of computing 
an element of the C(^) matrix, with all GPU cores 
running in parallel. The search of the minimum of 
C(^) is also performed by the GPU, thus reducing 
the transfers between GPU to CPU to the values 
of and 
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2.4-5. Monte-Carlo simulations 

Restricting and to belong to one of the 
DEM classes described in section [^XTl P{^^£,^) 
and P{£,^\£,^) are evaluated from Monte-Carlo 
simulations. For every combination of the two or 
three parameters defining (the ranges and res- 
olutions being given in 12.4.11) , 5000 independent 
realizations of the random and systematic errors 
are obtained. For each of the corresponding sets of 
six simulated AIA intensities, the inversion code 
returns the values of the parameters defining 
(equation ([5])) corresponding to the absolute min- 
imum of the criterion (equation (|15|) '). From the 
resulting 5000 we estimate the conditional prob- 
ability P{£,^\£,^) with a resolution defined by the 
sampling of the parameters. Integration over 
gives P{£,^) and using Bayes' theorem we obtain 
Pi^ie). 

3. Results: isothermal solution to isother- 
mal plasma 

In order to understand the fundamental prop- 
erties of the DEM inversion of the AIA data, 
we first applied the method to investigate the 
behaviour of the isothermal solutions to simula- 
tions of isothermal plasmas. The electron tem- 
peratures and emission measures of the plasmas 
are noted and EAl^ respectively. The cor- 
responding inverted quantities are noted T/ and 
EAf. The probabilities P{T^ , EM'\Tf , EM^) 
and P[T^ ,EMP\TI,EM^) are stored in matri- 
ces of dimension 4. To maximize the clarity of 
the results, and since the thermal content of the 
plasma is the main object of DEM analysis, we 
reduce the number of dimensions by fixing the 
emission measure of the simulated plasmas to be 
EM^ = 2 X 10^^ cm~^. Furthermore, the prob- 
abilities are always presented whatever the emis- 
sion measure by integrating them over EM^ , even 
though EM^ is of course solved for in the inversion 
process. 

The chosen E M^ is typical of non flaring ac- 



tive regions (e.g. IWarren et al. Il201ll) . Figure [3] 



shows as a function of and EM^ the num- 
ber of AIA bands in which a plasma produces 
more than 1 DN (detection threshold) and less 
than 11000 DN (saturation). The left panel is 
for isothermal plasmas, the right panel for Gaus- 
sian DEMs with (T^ = 0.51og(re). At the chosen 



EM^, and since we did not implement the detec- 
tor saturation in our simulations, we always have 
exploitable signal in all six AIA coronal bands, 
except below a few 10^ K. Conversely, solar struc- 
tures outside the white areas produce signal only 
in some of the six bands, unless spatial or tem- 
poral summation is used. Therefore, the results 
presented in the following sections correspond to 
optimum conditions outside of which the combina- 
tion of higher noise and possible lower number of 
valid bands will always lead to weaker constraints 
on the DEM. 

3.1. Three bands: EIT, TRACE, or low 
emission measures 

We first present inversion results using only 
three bands as an illustration of the situation en- 
countered with previous EUV imaging telescopes 
like EIT, TRACE or EUVI. The 17.1 nm and 
19.5 nm coronal passbands of EIT and TRACE 
have direct equivalents in AIA, but the Fe XV 
28.4 nm band does not. After comparison of 
its isothermal response (se e, e.g.. Figure 9 of 
iDelaboudiniere et al. 1995 ) with those of AIA 
(Figure [8]), we chose the Fe XIV 21.1 nm band 
as its closest AIA counterpart. The three bands 
configuration is also similar to having six bands 
and a low emission measure plasmgQ. Indeed, at 
5 X 10^^ m~^ and 1.5 x 10^ K, values typical of 
coronal loops, only three of the six AIA coronal 
bands produce more than 1 DN (see Figure [3]), 
the others providing only upper limit constraints 
to the DEM. 

Panel (a) of Figure [3] shows a map of the prob- 
abihtjH P{T^\T^). It is worth noting that, as ex- 
plained in section [O P(T/) and thus P{TP\T^) 
could be evaluated only because the limitation 
to simple parameterized DEM forms allowed the 
computation of P(r/|Tf ). The plot of P(r/) 
(and thus the horizontal structures in panel (a)) 
shows that some temperature solutions T/ are 
more probable than others for any plasma tem- 
perature Tf. In the case of real observations, 
this can be misinterpreted as the ubiquitous pres- 
ence of plasma at the most likely temperatures. 



■^For completeness, the plots for all combination of 
three to six bands are available on line at ftp.ias.u- 
psud.fr/cguennou/DEM_AIA_inversion/ 

^Defined as the probability for the solutions to lie between 
log Tc and log Tc + A log Tc . 
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This c aveat was already analysed by I Weber et al 
(|2005l ) in the case of the 19.5 to 17.3 nm TRACE 
band ratio and we will discuss it further in Paper 
II for multithermal plasmas. 

Both probability maps exhibit a diagonal from 
which several branches bifurcate. Below 2 x 10^ K 
and above 10^ K the diagonal disappears because 
since the bands have little sensitivity in these re- 
gions, the signal is dominated by noises and the 
inversion output is thus independent from the tem- 
perature. The general symmetry with respect to 
the diagonal reflects the equality P{T^\Tl) = 
P{T^\Tf)/P{T^). The diagonal is formed by the 
solutions T/ that are close to the input Tf , while 
the branches correspond to significant deviations 
from the input. In P{Tl\T^), these branches im- 
ply that two or more solutions T/ can be found 
for a same plasma temperature . Conversely, 
reading horizontally the P{T^\T^) image, a given 
temperature solution can be coherent with two 
or more plasma temperatures . The (b) and (c) 
plots give the probability of the solutions for two 
plasma temperatures. At = 3 x 10^ K, the so- 
lution may be T/ = 3 x 10^ K or 1.2 x 10^ K. At the 
typical coronal temperature = 1.5 x 10^ K, the 
inversion can yield 1.5 x 10^ K but also 2 x 10^ K 
or 10^ K. It is thus possible to incorrectly conclude 
to the presence of cool or hot coronal plasma while 
observing an average million degree corona. This 
ambiguity has far reaching implications since the 
detection of hot plasma i s one of the possible signa- 
tures of nano-flares (e.g. ICargill 1 ll994t iKlimchuk 
20061 ). Since by definition they correspond to the 



absolute minimum of the criterion, all solutions 
are fully consistent with the data given the uncer- 
tainties. One or more of the multiple solutions can 
be rejected only based on additional independent a 
priori information. For example, the high temper- 
ature solution corresponds to an emission measure 
of 4 X lO'^^ cm~^ (right panel of Figure[5]), which is 
extremely high considering the present knowledge 
of the corona. If no such information is available 
however, both low and high temperature solutions 
can still be correctly interpreted as also compati- 
ble with a 1.5 X 10^ K plasma with the aid of the 
P{TP\3 X IQS) and P(r/'|1G^) probability profiles 
(f) and (g). 

The reason for the formation of these branches 
is illustrated by Figure [H On both panels, the 
background image is the value of the criterion C((^) 



for a = 1.5 X 10^ K plasma as a function of 
Tc and EM. The absolute minimum of the crite- 
rion, the arguments of which are the inverted pa- 
rameters and EM^ (equation ([5])), corresponds 
to the darkest shade of grey and is marked by a 
white plus sign. The criterion is the sum of three 
components, one per waveband (equations (fT5|) 
and (fT7|) ). The three superimposed curves are the 
loci emission measure curves for each band b, i.e. 
the location of the {Tc, EM) pairs for which the 
theoretical intensities /^'' equal the measured ones 
I^^" ■ Below the loci curves, the criterion is almost 
flat because at lower emission measures the II^ 
are much smaller than the constant I^^" ■ Con- 
versely, the criterion is dominated by the Ith at 
high emission measures. The darkest shades of 
gray and thus the minimum of the criterion are 
located between these two regions. The two pan- 
els correspond to two independent realizations of 
the random and systematic errors. For each draw, 
the loci curves are randomly shifted along the EM 
axis around their average position. In the ab- 
sence of errors, the three loci curves would cross 
in a single point at the plasma temperature , 
giving a criterion strictly equal to zero. In the 
left panel, with random and systematic errors in- 
cluded, they do not intersect at a single point but 
the non-zero absolute minimum of C, where they 
are the closest together, is around . However, 
the criterion has two other local minima, around 
2 X 10^ K and around 10^ K, where two or three 
of the loci curves also bundle up. In the right 
panel, a different random draw shifts the curves 
closest together around the high temperature lo- 
cal minimum that thus becomes the new absolute 
minimum. For this 1.5 x 10^ K plasma, the inver- 
sion thus yields solutions randomly located around 
the several local minima with respective probabil- 
ities given by the profile of Figure IH^c). When 
scanning the plasma temperatures, the positions 
of the minima vary, thus building the branches in 
the probability maps. In addition, depending on 
their location the minima can be more or less ex- 
tended along one or the other axes, which results 
in a varying dispersion around the most probable 
solutions. 

Systematic errors are simulated with random 
variables while they are in fact identical for all 
measurements. Thus, the computed P(r/|T/') 
does not give the probability of solutions T/ for the 
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practical estimates of the calibration and atomic 
physics. In reality the output of the inversion 
is biased towards one or the other of the multi- 
ple solutions, but we do not know whether the 
calibration and atomic physics are under or over- 
estimated. Therefore, in order to deduce the prob- 
ability that the plasma has a temperature 
from an inverted temperature T/, we must ac- 
count for the probabilities of the systematics as 
defined in section [2.4.31 The randomization sam- 
ples their distribution, which ensures that the esti- 
mated P{T^\TI.) are the probabilities relevant to 
interpret T/. 

3.2. Six bands: AIA 

Figure [6] is the same as Figure |4j but now in- 
cluding the six AIA coronal bands in the analysis. 
Some secondary solutions persist at low probabili- 
ties but compared to the three bands case, most of 
the solutions are now concentrated on the diago- 
nal. This illustrates that the robustness of the in- 
version process increases with the number of bands 
or spectral lines. Comparison with Figure |4] quan- 
tifies the improvement brought by AIA over pre- 
vious instruments. Neglecting the low probability 
solutions, if independent a priori knowledge justi- 
fies the isothermal hypothesis, the six AIA bands 
thus provide an unambiguous determination of the 
plasma temperature. The temperature resolution 
of the inversion can be estimated from the width 
at half maximum of the diagonal. It varies over the 
temperature range between 0.03 and 0.11 \ogT^ . 
It is of course be modified if we assumed differ- 
ent uncertainties on the calibration and atomic 
physics than the ones chosen in section [2.4.31 We 
tested the sensitivity of the temperature resolu- 
tion to the level of uncertainties cr^, from 10 to 
55%. The higher the uncertainties, the lower the 
temperature resolution and the more probable the 
secondary solution^ For an estimated temper- 
ature of 1 MK, the temperature resolution of 
the inversion varies between 0.02 for 10% error and 
0.08 logTjf for 55% error. In the worst case, for 
55% errors, the temperature resolution decreases 
to 0.2 \ogT^ for the temperature interval between 
0.5 and 0.9 MK. At 1 MK, the resolution is pro- 
portional to the uncertainty level with a coefficient 



of 0.15 (AT/" - 0.15 <). 

Since by definition our method always finds 
the absolute minimum of the least square crite- 
rion of equation (jlSp . the derived temperature 
resolution is an intrinsic property of the data 
and not of the inversion scheme. It is the re- 
sult of the combination of the random and sys- 
tematic errors and the shapes the contribution 
functions. Its value is direct l y com parable with 
the findings of iLandi et al. I (j201iri . These au- 



thors showed that the tempe rature re s olutio n of 
the MCMC code of .Kashvap fc Drake I (|l998[ ) ap- 
plied to isothermal plasmas is 0.05 log T. Their 
tests were made on simulated observations of a 
10^ K plasma in 45 isolated spectral lines with 
20% random errors. Assuming that the MCMC 
method does converge towards solutions consistent 
with the limitations of the data, the fact that the 
temperature resolution is comparable for 6 AIA 
bands and 45 spectral lines suggests that, in the 
isothermal case, it is driven by the uncertainties 
level rather than the number of observables. This 
conclusion is consist ent with the isotherm al limit 
of Figure 6 of Tandi fc KUmchuk I (|2010l ). 



"The corresponding probability maps are available on line 
at ftp.ias.u-psud.fr/cgucnnou/DEM_AIA_invcrsion/ 



3.3. Residuals and goodness of fit test 

The probability maps presented in the above 
sections are valid for a given hypothesis on the 
plasma DEM distribution, but they would be use- 
less without a test of its validity. The pertinence 
of the DEM model chosen to interpret the obser- 
vations can be assessed by analyzing the distri- 
bution of the sum of squared residuals defined by 
equations and If applying our inversion 
scheme to real data, we could compare the result- 
ing residuals to the distribution derived from sim- 
ulations for a given DEM model and thus quantify 
the probability that the data is consistent with the 
working hypothesis (e.g. isothermal or Gaussian). 

The solid line histogram of Figure [7] shows 
the distribution of values corresponding to the 
plots of Figure |6l The distribution is close to a de- 
gree 4 distribution (solid curve) although not a 
perfect match, with a peak shifted to the left and 
an enhanced wing. The most probable value of the 
squared residuals is ^ 1.73 and 95% of them are 
comprised between and 15. Whatever the actual 
plasma DEM, any inversion made with the isother- 
mal hypothesis and yielding a value in this 
range can thus be considered consistent with an 
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isothermal plasma given the uncertainties. This 
iso thermality test is similar to that recommended 



bv lLandifcKlimchukI ^oE): identifying our 
to their i^min and our maximum acceptable 
to their AF. This does not imply however that 
isothermality is the only nor the best possible in- 
terpretation of the data, for different DEMs can 
produce similar values. The discrimination be- 
tween DEM models will be discussed in Paper II. 

The properties of the empirical distribution 
of squared residuals can be explained as fol- 
lows. Since we simulated observations of a purely 
isothermal plasma, an isothermal model can al- 
ways represent the data. Without errors, there 
would always be one unique couple (T/, EM^) 
corresponding to six intensities perfectly matching 
the six AIA observations, thus giving zero resid- 
uals. With errors, if we forced the solution (T/, 
EM^) to be the input {T^ , EM^), the summed 
squared residuals resulting from a number of ran- 
dom draws should have by definition the Proba- 
bility Density Function (PDF) of a degree six 
distribution (dotted curve of Figure [7]), for we have 
six independent values of 7^'"' — 7^'* and we nor- 
malized the residuals to the standard deviation 
(7^ of the uncertainties. But since we solve for 
two parameters (T/, EAI^) by performing a least 
squares minimization at each realization of the 
errors, the solution is not exactly the input {Tf, 
EM^) and we should expect a PDF with two de- 
grees of freedom less (dashed curve). Instead of 
being a pure degree 4, the observed distribution 
is slightly shifted toward a degree 3 because of 
two factors. First, the errors are a combination 
of Poisson photon noise and Gaussian read noise, 
while the x? distribution is defined for standard 
normal random variables. Second, as discussed 
below, the six residuals are not completely inde- 
pendent. 

Figure [5] shows the response functions R}j{Tf,) 
of the AIA bands to isothermal plasmas with elec- 
tron temperatures from 10^ to 3 x 10^ K for a 
constant electron number density of 10^ cm^'^. 
For each band, the thick curve is the total re- 
sponse, and the labeled thin curves are the partial 
responses for the ions that contribute the most 
for at least one temperature. The fraction of the 
total response not accounted for by those domi- 
nant ions is shown below each main plot. Ion- 
ization stages common to several bands are found 



across the whole range of temperatures. O V dom- 
inates the response at 2.5 x 10^ K in the 17.1 nm, 
19.3 nm and 21.1 nm bands. Around 1 MK, Fe IX 
is found in the 17.1 nm, 19.3 nm and 21.1 nm 
responses, and Fe X contributes to the 94 nm, 
21.1 nm and 33.5 nm bands. At 2 MK, Fe XIV 
is common to the 21.1 nm and 33.5 nm bands. 
This is c onsistent with t h e ana lysis of the AIA 
bands bv lO'Dwver et al. I (|20inl ). Because of this 
redundancy, the response functions tend to have 
similar shapes in the regions of overlap, resulting 
in a correlation between the residuals. 

4. Summary and conclusions 

By restricting the solutions to functional forms 
described by a limited number of parameters, we 
obtained a complete statistical characterization of 
the DEM inversion. Even though they are not ex- 
pected to accurately describe real coronal proper- 
ties, these simple DEM distributions can nonethe- 
less model a wide range of plasma conditions. The 
results presented in this series of papers can thus 
be fruitfully used to demonstrate many important 
properties and guide the interpretation of the out- 
put of generic DEM inversion codes. We illus- 
trated the method by applying it to the six coro- 
nal bands of the AIA telescope. In this first pa- 
per, we limited ourselves to isothermal plasmas 
and isothermal solutions. 



The case presented in section [3.11 demonstrates 
the existence of multiple solutions if the number 
of bands is limited either by design of the instru- 
ment or by lack of signal. However, since our 
method provides the respective probabilities of the 
multiple solutions, it is possible to properly in- 
terpret the solutions as compatible with several 
plasma temperatures. Even if some of these prop- 
erties have been illustrated in case studies, we pro- 
vide here a systematic analysis of a wide range of 
plasma parameters. The computed distribution of 
squared residuals can be used to test the coherence 
of real AIA data with the isothermal hypothesis. 
This type of analysis can also be help to deter- 
mine the optimum data acquisition parameters for 
AIA (e.g. spatial binning and exposure time) en- 
suring that no secondary solution is present. In 
section 13.21 we showed that, with enough signal, 
the six AIA coronal bands provide a robust recon- 
struction of isothermal plasmas with a tempera- 
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ture resolution comprised between 0.03 and 0.11 
logTg. The comparison with the three bands case 
gives a quantification of the improvement brought 
by the new generation of instruments. The same 
method can be appHed to other instruments with 
different response functions and different numbers 
of bands or spectral lines. This naturally requires 
the computation of the corresponding probability 
matrices and distribution of residuals. 

The temperature resolution, and more gener- 
ally the details of the probability matrices pre- 
sented in sections 13.11 and 13.21 depend on the am- 
plitude and distribution of the random and sys- 
tematic errors. We found the resolution to be 
proportional to the uncertainty level (at 1 MK, 
at/' ^ 0.15 cr^). We simulated plasmas with 
high emission measures typical of active regions. 
Depending on the temperature, either the photon 
noise or the uncertainties on the calibration and 
atomic physics dominate. The illustrated prop- 
erties of the inversion, from the multiplicity of 
solutions to the temperature resolution, are thus 
driven by both random systematic errors. While 
the photon noise can be reduced by increasing the 
exposure time or binning the data, reducing the 
systematics requires better atomic data and pho- 
tometric calibration, which is not trivial. 
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Fig. 1. — Principle of the method used. Reference theoretical intensities are tabulated using CHIANTI 
for different parameterized DEM functional forms (Dirac, Gaussian, top hat). A random variable is added 
to represent the uncertainties on the calibration and atomic physics. For a given plasma DEM AIA 
observations are simulated in a similar way. A criterion is minimized to find the DEM that best 
matches the simulated observations. By scanning the parameters defining the probabilities ^'(^^l^^) 
and P(C^IC^) are built from a large number of draws of the random variables. These probabilities and the 
corresponding distributions of values give a complete characterization of the inversion for the chosen 
DEM forms. 
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Fig. 2. — Normalized maximum of the response functions Rb(Te,ne) of the six AIA coronal bands as a 
function of electron number density. Only the 9.4 nm band is independent on the density as assumed in the 
DEM analysis. The other functions vary by up to 35% in the range of densities plausible in the AIA field of 
view (dashed vertical lines). This effect induces systematic errors in the DEM inversions. 
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Central temoerature (K) 

Fig. 3. — Number of AIA bands in which the signal is comprised between 1 DN (detection threshold) and 
11000 DN (saturation) as a function of temperature and emission measure, for standard exposure times. 
Left: isothermal plasmas ; right; Gaussian DEMs with a = 0.5 log(Te). Only solar structures falling in the 
white regions produce exploitable signal in all six AIA coronal bands. The regions corresponding to five valid 
bands are labelled with the wavelength of the missing one. If exposure times were increased, the boundaries 
of all regions would be shifted towards smaller omission measures. If several images were summed up to 
overcome saturation, the upper boundaries would be moved upwards. 
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Fig. 4. — Probabilities of the isothermal solutions for observations of an isothermal plasma with three of the 
AIA coronal bands (17.1, 19.3 and 21.1 nm). (a): reading vertically, conditional probability P{T^\T^) that 
the inversion yields T/ for a given plasma temperature . (e): reading horizontally, probability P{T^\Tj.) 
that the plasma has a temperature for an inverted value T/. (e) is obtained by normalizing (a) to (d) 
the unconditional probability P(T/) that the inversion yields T/ whatever the plasma temperature, which 
is obtained by integrating (a) over . The branches bifurcating from the diagonal reveal the existence of 
multiple solutions. The probability profiles (b) and (c) show for example that 3 x 10^ K or 1.5 x 10® K 
plasmas can be measured at 3 x 10^ K, 1.5 x 10® K or 10^ K. Vice versa, the profiles (f) and (g) can be used 
to properly interpret 3 x 10^ K and 10^ K inversions as also both compatible with a 1.5 x 10® K plasma. See 
section [XT] for details. 
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Fig. 5. — The least squares isothermal criterion (Equation (fTSj)) for a simulated isothermal plasma at 
= 1.5 X 10^ K. The two panels correspond to two of the random draws used to build Figure 01 The 
loci curves for the three components of the criterion are superimposed, and the white plus signs mark the 
location of its absolute minimum. Both solutions are fully consistent with the simulated data given the 
uncertainties. 
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Fig. 6. — Same as Figure 2] for the six AIA coronal bands. The determination of the temperature of the 
simulated isothermal plasma is now unambiguous. From the width of the diagonal we deduce the resolution 
of the isothermal inversion to be about 0.05 log(Tc). 
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Fig. 7. — The observed distribution of the sum 
of the squared residuals (solid histogram) differs 
somewhat from the expected degree 4 dis- 
tribution (solid curve). It is slightly shifted toward 
a degree 3 (dashed curve), which can be explained 
by a small correlation between the six AIA coronal 
bands. 
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Fig. 8. — Isothermal response of the six AIA coronal bands between 10^ and 3 x 10^ K. For each band, the 
thick curve is the total response and the thin curves are the partial responses for the ions that contribute 
the most at at least one temperature. The fraction S2 the total response accounted for by those dominant 
ions is shown below each main plot. Computations for an electron number density of 10^ cm~^. 



